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Abstract. The problem of the interplay between normal and anomalous scaling in 
turbulent systems stirred by a random forcing with a power law spectrum is addressed. 
We consider both linear and nonlinear systems. As for the linear case, we study passive 
scalars advected by a 2d velocity field in the inverse cascade regime. For the nonlinear 
case, we review a recent investigation of 3d Navier-Stokes turbulence, and we present 
new quantitative results for shell models of turbulence. We show that to get firm 
statements is necessary to reach considerably high resolutions due to the presence of 
unavoidable subleading terms affecting all correlation functions. All findings support 
universality of anomalous scaling for the small scale fluctuations. 
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1. Introduction 

The understanding of the small scale statistics of turbulent systems is a problem of 
considerable interest [I]. With turbulent systems, we mean both the dynamics of 
velocity fields in high-Reynolds number flows, and the advection of scalar or vector 
fields by turbulent flows as, e.g., temperature or magnetic fields. In the last few years, 
much progress has been reached both in experimental O El 13 E] and numerical [HI E] 
investigations of turbulent systems. We have now plenty of observations showing that 
two properties generally holds. First, turbulent systems are intermittent El El , 
as quantified by the anomalous scaling of moments of field increments. Second, 
the anomalous scaling exponents display universality with respect to the boundary 
conditions and to the large scale forcing mechanisms [3 El- Some subtle points arise 
when isotropy is broken at large scale by the injection mechanisms |l() | IHllT H ll2j. In that 
case, for universality to hold it is required that anisotropic contributions are subleading 
with respect to the isotropic one. However, evidences for universality of both isotropic 
and anisotropic scaling exponents have been found [T^ I14j. 

Nonetheless, the mechanisms responsible for anomalous scaling and universality 
in turbulent systems are still not fully comprehended with the remarkable exception 
of linear problems such as the passive transport of a field, for which intermittency and 
universality have been systematically understood (see Ref. for an exhaustive review). 
In this context, in particular for the class of Kraichnan models |TS1, closed equations for 
the correlation functions can be derived. These are linear partial differential equations, 
whose homogeneous solutions (zero modes) generally exhibit anomalous scaling. On the 
other hand, the inhomogeneous solutions, constrained by the external forcing, possess 
dimensional (non-anomalous) scaling. Universality results then from the decoupling 
between the zero modes scaling and the forcing. Remarkably, the zero modes can be 
interpreted as statistically preserved structures, i.e. functions that do not change in time 
once averaged over the velocity field realisations and particle trajectories. This allowed 
for successfully testing the entire picture also for the advection by realistic velocity fields 
[T71ITH1, and in the context of shell models for passive transport P^ . 

On the theoretical side, it is tempting to export the concepts issuing from linear 
turbulent systems to nonlinear ones, i.e. to see whether the same mechanisms for 
anomalous scaling and universality are at work also in nonlinear hydrodynamic systems 
such as Navier-Stokes turbulence. A possible test case to probe universality with respect 
to the forcing mechanisms is to study turbulent fluctuations stirred at all scales by a self- 
similar power-law random forcing. Indeed, the non- analytical properties of the forcing 
may alter the energy exchange between the turbulent field fluctuations. For instance, in 
some cases, the energy injection mechanism may prevail over the energy cascade process, 
modifying the inertial range physics. 

This problem was pioneered by means of Renormalisation Group (RG) methods 
[20, ^|22l, in the late '70s, for the d-dimensional Navier-Stokes equations. In the RG 
calculations, for d = 3, the forcing spectrum is chosen as Ef{k) ~ k^'^ with y playing 
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the role of a small parameter in perturbative expansions. Unfortunately the interesting 
physical case, obtained for y = 4 and corresponding to the Kolmogorov spectrum for the 
velocity field, lies in a range where convergence of the RG expansion is not granted |25j . 
Notwithstanding extensions of the RG formalism to ?/ ~ 0(1) values have been tempted 
by different approaches [231 El], the problem is still open. Recent numerical simulations 
tried to shed light on this issue [23 I2Z] but, because of the limited resolution, their 
results are not conclusive. 

The aim of this paper is to investigate the general issue of the small scale statistical 
properties of linear and nonlinear hydro-dynamical systems, in the presence of stirring 
acting at all scales with a power-law spectrum. A systematic study of the scaling 
behaviour at varying the forcing spectrum allows indeed for understanding the interplay 
of dimensional and anomalous scaling in turbulent fields. 

In Sec. 121 as an instance of linear problems, we consider passive scalars stirred 
at all scales by a power-law forcing, and advected by a 2d turbulent velocity field. 
Our results confirm the universality scenario originating from the zero modes picture 
which predicts two distinct regimes. Forcing- dominated regime: the scaling of low 
order structure functions is non-anomalous, with exponents dimensionally related to the 
forcing spectrum; for the higher order moments, scaling is anomalous and dominated 
by the zero modes. Forcing- subleading regime: the dimensional scaling related to the 
balance with the forcing is subleading, at any order, with respect to the anomalous one, 
similarly to the case of a standard large scale injection. Hence, anomalous scaling is 
observed for any order statistics. A subtle technical point revealed by the passive scalar 
case is the existence of many important power-law terms that contribute to the scaling 
properties. As clarified in Sec. [21 for both regimes, to disentangle the authentic scaling 
behaviours, it is necessary to take into account the leading as well as the subleading 
terms. Concerning nonlinear systems, in Sec. [21 we review the numerical study of 
the 3d Navier-Stokes equations done in Ref. ^2!, where due to natural limitation in 
the resolution only semi-quantitative results were obtained. Then, to overcome this 
difficulty, we consider shell models for turbulence. These are nonlinear models that 
maintain most of the richness of the original problem but allow for reaching higher 
Reynolds numbers. Shell model results coherently fit with the passive scalar and NS 
ones. Conclusions follow in Sec. HI 

2. Linear dynamics: passive scalar transport 

The evolution of a passive scalar field 6 advected by an incompressible fiow v is governed 
by the advection-diffusion equation : 

dtO + v-Ve = kAO + f , (1) 

where k is the molecular diffusivity. The standard phenomenology is as follows. Scalar 
fiuctuations are injected at large scale by the source term /; then, by a cascade 
mechanism induced by the advection term, they reach small scales where dissipation 
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takes place balancing the input, and holding the system in a statistically stationary 
state. 

Since the problem is passive, it can be studied also with a synthetic velocity 
field mimicking some of the features of realistic turbulent fiows. Much attention 
has been recently devoted to the Kraichnan model jTH], where i; is a Gaussian, 
incompressible, homogeneous, isotropic and 5-correlated in time field of zero mean. The 
only reminiscence of real turbulent fiows is the existence of a scaling behaviour that, 
due to the assumed self-similarity, is fixed by a unique exponent, i.e. {{6rV ■ f)^) ~ 
with < ^ < 2 (where SrV = v{x + r) — v{x)). For the sake of analytical control, the 
forcing / is also taken as a Gaussian, homogeneous, isotropic field of zero mean and 
with correlation function: (/(r, t)/(0, 0)) = 5(t)jF(r), where J-'{r) decays rapidly for 
r ^ Lf, identifying Lf as the large scale of the problem. 

Thanks to the above simplifying assumptions, a closed equation for the generic p- 
point correlator, Cp{r,t) = {9{ri,t)9{r2,t)...9{rp,t)) , can be derived. It formally reads 

dtCp = -MpCp + J^® Cp.2 , (2) 

where T is the forcing spatial correlation defined above, and M.p is a linear differential 
operator coming from the diffusion and the advection terms (see Ref. ^H] for details). 
The stationary solution of Q satisfies the equation M.pCp = jF(g)Cp_2. As usual, a linear 
equation is solved by the superposition of the solution of the associated homogeneous 
equation, J^pZp = 0, and the solution of the inhomogeneous one, which will be denoted 
as Cp. In the inertial range of scales both Cp and Zp display a scaling behaviour, i.e. 

Cj(Ar) ~ A^"'"Cj(r) and Zp{Xr) X^' Zp{r) , with r = (n, . . . , r^) , (3) 

where the scaling exponents Cp*™" = p(2 — ^)/2 are fixed by dimensionally matching 
the inertial operator with the forcing, while the (p, not constrained by any dimensional 
requirements, are independent of the forcing. One is usually interested in the p-point 
irreducible component of the correlation function, Cp, that is the structure function 
Sp{r) = {{6r0y). Still in the framework of Kraichnan model, the zero modes scaling 
exponents contributing to Sp have been calculated perturbatively in [2H1 120] • It has been 
found that (p < Cp*™ for p > 2, i.e. the leading contribution to the structure function 
scaling is anomalous. 

The zero modes dominance scenario explains both the anomalous scaling exponents 
and their universality. Indeed the forcing does not enter the definition of Zp, but fixes 
only the multiplicative constants in front of Cp and Zp, necessary to match the large scale 
boundary conditions. This implies that, even though scaling exponents are universal, 
probability density functions (PDF) of scalar increments are not. The above picture 
based on zero modes is now recognised to apply in all linear hydrodynamic problems. 
As confirmed by investigations of passive advection by realistic velocity fields [T^ ITTj. 
and of shell models for passive transport [T^ . 

Let us now focus on the main issue of this work: the scaling behaviour in the 
presence of forcing fiuctuations directly injected in the inertial range of scales. We 
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consider the case of a source term f{r,t) which is a random Gaussian scalar field, with 
zero mean, white-in-time and characterised by the spectrum 



In the following, it is always assumed the presence of ultraviolet and infrared cutoffs for 
the forcing, i.e. expression holds only in the range ki < \k\ < k2 oi wave-numbers, 
where ki and ^2 are of the order of the inverse of the largest scale in the system and of 
the dissipative scale, respectively. 

Remaining in the framework of the Kraichnan model, Eq. (0) still holds, and its 
stationary solution is again the superposition of the homogeneous solution, 2p, which 
remains unchanged {viz. with the same scaling exponents) and of the inhomogeneous 
one, Cp , whose scaling exponents now depend on the slope of the forcing spectrum, 13. 
Indeed by dimensional reasoning, we have Cp^'^iP) = p(2 — C, — /5)/2. 

Two regimes can be identified. If /3 < 0, the scaling of the inhomogeneous solution is 
always subleading with respect to the anomalous one. The scaling exponents measured 
from the structure functions are the same as those obtained with the standard large 
scale forcing. On the other hand, for (3 > 0, there exists a critical order pc such that 
for p < Pc- Cf'^i.P) < Cp) and for p > p^. Cp < Cp*'"(/^)- other words, the scaling 
behaviour of the low-order structure functions is non-anomalous and dominated by the 
forcing. The appearing of anomalous scaling for p > Pc can be understood due the fact 
that (p as a function of p is concave. Moreover, it is known numerically fH] and to 
some degree analytically [HHI that in the Kraichnan model the scaling exponents above 
a certain order saturate to a finite value, Coo? which makes even more intuitive the 
existence of a finite pc- Some subtle points may arise for positive /3 larger than 1. In 
such a case, the strong UV components of the forcing spectrum may allow a matching 
with zero modes, disregarded for a large scale forcing [221, exploding at small scales. 
Finally, the case (3 = is marginal, because the exponents Cp*"* coincide with the large 
scale prediction, up to possible logarithmic corrections. 

Checking the validity of the above predictions in numerical simulations of a realistic 
fiow is interesting for two reasons. First, it is a further demonstration, in the Eulerian 
framework, of the zero modes picture for anomalous scaling and universality in linear 
problems, which was previously assayed with Lagrangian studies ^Tj. Second, it offers a 
controlled testing ground for interpreting some aspects of the nonlinear hydrodynamics 
which will be considered in the next section. 

As an example of realistic velocity field, we consider two-dimensional incompressible 
Navier-Stokes equations in the inverse cascade regime [SHI HI IHl] : 



The terms in the r.h.s. of Eq. (0) have the following meaning: P is the pressure field; /„ 
injects kinetic energy at scale Lq] uAv (where u is the viscosity) dissipates enstrophy 
at small scale; —av removes energy at the large scales allowing for a statistically 
stationary state. In the inverse energy cascade regime, the velocity statistics is self- 
similar (not-intermittent), with {{5rV ■ rY) ~ r^/^^ but temporal correlations are non 



Ef{k) = T,k{\f{k)\')^k-~'^^. 



(4) 



dtv + V ■ V V = — VP + TjAv — av + 



(5) 
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Figure 1. Instantaneous snapshots of the scalar field 9 for runs (a), (b), and (c); 
intensities are coded in grey-scale. Note that (a) and (b) are qualitatively very similar, 
while (c) is characterised by richer fluctuations at the small scales. The three snapshots 
are obtained by integrating Eqs. and Q by means of a standard 2/3-dealiased 
pseudospectral code in a doubly periodic square domain 2tt x 2tt with 1024^ grid points. 
As for the velocity field, in Eq. Q the viscous term has been replaced by a hyper- 
viscous term of order 8, in order to force at the very small scales. The velocity forcing 
fy is chosen as a Gaussian, incompressible, homogeneous, isotropic and ^-correlated in 
time 2d field, of zero mean and concentrated around the small scale Lq (of the order 
of few grid points). The friction coefficient a is tuned in such a way that energy is 
removed at a scale rjfr ~ (tj'^a^'^l'^ of the order of the box size. In Eq. to have 
a large inertial range, the diffusive term has been replaced by a bilaplacian with k 
tuned to have the dissipative scale > Lq. As for the scalar, in run (a) we used a 
white- in-time, random Gaussian forcing concentrated at scale i/ of the order of the 
box size. In the runs (b) and (c) we used the forcing Q by choosing ki — 2 and ^2 in 
such a way that k2 < 2'iT/rd. For the sake of comparison, we imposed in all runs the 
same scalar energy input. For each run, we collected 80 frames separated by about one 
large scale eddy turnover time measured as rjfr/ \J (u^). See [3^] for more details on 
the numerical procedure. 

trivial. Moreover, precise numerical measurements of passive scalars with large scale 
forcing ^SjEI have shown that: ^2 = 2/3 while for p > 2 the exponents are anomalous, 
and saturation is observed for p > 10 with (oo — 1-4. It is noteworthy that saturation 
seems to be generic in passive scalars, as found also in experiments |3IJ. From a physical 
point of view, it means that the strongest scalar fluctuations are statistically dominated 
by front-like structures, i.e. huge variations of the field 6 at very small scales. 

Once a power law forcing is considered, the predicted dimensional scaling 
is Cp''"^{P) = p(2/3 — P)/2. Therefore, the two above described regimes -statistics 
dominated by zero modes or by the forcing- should appear for /? < and P > 0, 
respectively. We performed three sets of direct numerical simulations (DNS) of 
Eqs. (H)),©: for run (a), we used a standard large scale forcing; for runs (b) and (c), we 
considered a forcing as in Eq. (jlj) with P = —0.3 and P = 0.3, respectively. More details 
on the DNS are given in the caption of Fig. ^ where we show the snapshots of the scalar 
field 9 for all runs. Already at a first sight, it is possible to observe that runs (a) and 
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Figure 2. Log-lin plot of local slopes of the second order structure function for the 
three runs. The dotted straight line indicates the value 2/3. Note that in run (a) the 
slope are constant for about one decade. On the other hand, the slopes for run (b) 
and (c) are non constant in the entire range. It is worth remarking that for power-law 
forcing there are two sources of finite size effects. The first is the presence of two 
power-laws, see text. The second is that the scaling in real space of the forcing two- 
point correlation is strongly affected by corrections (induced by the Bessel function) 
due to fact that we generate the forcing in Fourier space. 



(b), where the zero modes dominate over the inhomogeneous solutions at any order, 
display qualitatively similar features: the field is organised into large scale structures 
or plateaux characterised by a good mixing, separated by sharp fronts. Differently in 
run (c), which corresponds to the case of forcing dominated statistics (with Pc>6), the 
scalar fluctuations develop a wider range of scales, and the most evident large scales 
plateaux disappear. 

In order to make the above observations more quantitative, we studied the scaling 
of the structure functions S2p{r) = {{6r9Y^) (odd orders are zero due to the isotropy of 
the forcing). For run (a), we found a very good scaling range of about one decade, which 
allowed us for accurately measuring (p up to p ~ 10 — 12. The quality of the local slopes 
and the measured values agree with previous investigations jlHI, which were performed 
with larger statistics and higher resolution. For the power-law cases, runs (b) and (c), 
we could not flnd evidence of a good scaling as in run (a). This is observed already 
from the second order structure function 5*2, whose local slope is plotted in Fig. |21for 
the three runs. 

We understand the poor scaling behaviour as due to the competition between the 
scaling of the anomalous part and the forcing dominated one. This becomes clear by 
looking at the scalar flux, Ilg{r) = {[6rV ■ f](5r^^)^)- For this quantity under rather 
general hypothesis, an exact analytical prediction can be derived for a generic forcing, 
through the Yaglom equation For a power law forcing / as in (0]), we obtain: 

Ueir) = {[5rv ■ f] (5,0)2) ~ cir + car^-^ , (6) 
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where we have kept only the first two leading contributions. It is important to notice 
that constants Ci and C2, which depend on the space dimension and on the details of the 
forcing spectrum, turn out to have opposite signs. In (jH]), in addition to the standard 
linear term, there is the first leading term induced by the forcing (jD). As one can see, 
the critical value (3 = naturally arises in the flux expression to separate the inertial 
dominated {(3 < 0) from the forcing dominated (/5 > 0) regime. Indeed, for j3 < the 
spectral flux is dominated by the low wavenumber components of the forcing spectrum 
and saturates, in Fourier space, to a constant value as a function of k. Scalar fluctuations 
are transferred down-scale via an intermittent cascade and the statistics is dominated by 
the anomalous scaling, the forcing contribution even if present at all scales is subleading. 
For (3 > 0, the scalar spectral flux no longer saturates to a constant: the direct input of 
energy from the forcing mechanism affects inertial range statistics in a self-similar way, 
down to the smallest scales where dissipative terms start to be important. 

In addition to Il0{r), we also studied the third moment of the flux variable 
Ilg{r) = {[{SrV ■ r){6r9)^]^) whose scahng is anomalous for the large-scale forcing case, 
i.e. ng(r) ~ r^ "^ (see also [HI)- Results are summarised in Fig. El 




Figure 3. Left: structure functions of passive scalar flux variables. Blue crosses: run 
(a); green squares: run (b); red circles: run (c). Fits are made in the inertial range with 
the expression given by the Yaglom equation ©. Insets: same analysis is performed 
for compensated ratios in the case of run (b) (top inset) and run (c) (bottom inset). 
We show the best compensation obtained with both power laws (same symbols as in 
the body of the figure) in Eq. ©, and with only the leading term (dotted lines). Right: 
the same of left panel for the third order moment of the flux variable. For the large 
scale forcing we use the single power-law, r^ "*, to fit the data, while for the other two 
cases we used a power law superposition of the two terms of Eq. Q . For all fits with 
a power-law superposition, the two terms turns out to have opposite signs. Insets: 
compensations with a single (leading) or two power laws for run (b) (top inset) and 
run (c) (bottom inset). Here also the best compensation is obtained with two power 
laws. 

Concerning the flux Ile{r), for the large scale case run (a), the scaling is good and 
the Yaglom relation, Ile{r) ^ —2eer (being eg the scalar dissipation), holds for about 
one decade (Fig. 01 left panel). On the other hand, for the power-law forced cases, the 
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scaling is poorer, making the identification of tlie exponents less trivial. Indeed, a scaling 
behaviour can be properly identified only by taking into account both the leading and 
subleading terms in Eq. (jHl). The third order moment is affected by the same problem. 
So that, on the basis of relation obtained for Uq, we tried to fit it with the following 
scaling ansatz: 

n^(r) ~ Ar2-^ + 5r3(i-'^), (7) 

where the exponent 2.4 in the first term is associated to the anomalous contribution as 
measured from run (a). Notice that for run (b) the exponent induced by the forcing is 
subleading, i.e. 3(1 — /3) = 3.9, while for run (c) it becomes leading 3(1 — /3) = 2.1. 
From the insets of Fig. IHl it is clear that a well defined scaling behaviour is recovered 
only after having taken into account the two power-law contributions. In all cases the 
anomalous scaling exponents have values in agreement with those obtained for a large 
scale forcing, although the interplay of different power laws, when the the forcing is as 
(j3)), makes the identification of the exponents very difficult even at considerably high 
resolution. 




Figure 4. Top: PDF of scalar increments, 6rO, for runs (a), (b) and (c). For each run, 
we plot the PDF for three different separations r inside the inertial range. The curves 
have been normalised to have unit variance. The dotted line is a Gaussian distribution 
with unit variance for comparison. Bottom: the same but zooming in, to highlight the 
core behaviour. Note that for run (c) the PDF's core approaches a Gaussian with a 
weaker, if not absent, dependence of r. 

Due to these difficulties in measuring the exponents, we looked directly at the PDF 
of scalar increments P{SrO). In Fig. EJ we show plots of normalised PDFs of 6r0 for 
different choices of the scale r in the inertial range. For runs (a) and (b), the PDFs do 
not display any rescaling property at changing the scale r, neither in the core nor in the 
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tails, suggesting intermittent behaviour in the full statistics. Note also that the PDFs 
in the two runs are different (only scaling exponents are universal while PDFs are not). 
Conversely, for run (c), where the forcing is dominant, the PDFs display a fairly good 
rescaling in the core (which governs the low order statistics), while the tails still do not 
collapse. This agrees with the existence of a critical order, pc, above which anomalous 
scaling appears even in the forcing-dominated case. As discussed earlier, all differences 
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Figure 5. PDF of scalar increments rescaled to test Eq. (jSJ with Coo — 1-4 for all runs 
(a), (b) and (c). For plotting purposes, curves of runs (b) and (c) have been shifted 
by a factor 10'^ and 10^, respectively. 



associated with the two regimes should disappear in the high order statistics, i.e. for 
p > Pc when anomalous scaling exponents imposed by the zero modes should show up 
irrespectively of the forcing. In testing this point, the presence of saturation for large 
p comes into help, indeed it entails that, for large excursions, the PDFs must approach 
the form [TH] : 

P{S,e) = r^-Q (^] (8) 

\ "rms / "rms 

with 9rms = \[(^ s-iid Coo — 1-4. Such a prediction is tested and confirmed by the 
collapse in the PDF tails shown in Fig. El which tells us that Coo is asymptotically 
approached independently of the forcing. 

The above results provide support to the validity of the zero modes picture 
beyond the boundaries of the Kraichnan model and in agreement with the Lagrangian 
investigations [T7j . 



3. Nonlinear systems 

Observations from the linear problem strongly indicate that the notions of anomalous 
scaling and universality are closely linked. When we enter the nonlinear world, we do 
not have any longer a reference theory. Here, we shall rather formulate some working 
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hypothesis suggested by the results from the hnear systems and verify their consistency 
with two cases studied, namely: 3d Navier-Stokes turbulence, reviewing the work first 
presented in and shell models jHE] for turbulence. 

3.1. The 3d Navier-Stokes problem 

In the 3d case, we consider random forcing defined by a two-point correlation function 
that in Fourier space reads 

{Mk,t)f,{k',t'))cxk'-yp,,{k)S{k + k')s{t-t'), (9) 

where Pij{k) is the projector assuring incompressibility and the forcing spectrum 
behaves as Ef{k) ~ k^^^. Hereafter, we use the definition y = 4 — (3, referring to 
the classical notation of the problem as introduced in [201 Correspondingly, the 
critical value (3 = 0, separating the two regimes in the linear case, is now yc = 4. 

The relative importance of the stirring on the small scales can be varied by 
tuning the slope value y from y ~ (meaning strong input at all scales) to y — > cxo 
(corresponding to a quasi large scale forcing). Also in NS turbulence, a physical insight 
can be taken from the behaviour of the energy flux, which is constant through scales 
up to logarithmic corrections for the subleading forcing case, y > yc] while it becomes a 
scale dependent function for y < y^., overcoming the energy cascade mechanism. 

The case of strong input at all scales, y ~ 0, was originally investigated in Ref. [20] 
by means of a RG approach, leading to the following expression for the energy spectrum 
E{k) ~ k^~'^y^^, in the domain t] <^ k~^ ^ where r] and Lf are the viscous scale 
and large scale of the system, respectively. It is worth noticing that such a prediction, 
which results also from dimensional analysis leads to the Kolmogorov spectrum 
E(k) ~ k~^^^ for y = A, i.e. quite far from the perturbative region where the RG 
calculations are under control. Here we want to study the same problems addressed for 
the linear case, i.e. whether fluctuations are sensitive to the injection mechanism for 
any y and if one observes anomalous scaling for y < 4, when the forcing is dominant. 
Since, at least with the present knowledge, RG perturbative methods starting at ?/ ~ 
cannot control anomalous corrections, only numerical simulations at finite y values, can 
possibly give some answers. 

In Ref. jSHj a first numerical investigation of 3d incompressible Navier-Stokes 
problem was performed. The authors made a set of DNS, varying the spectrum slope 
from y = 3 to y = 8 ai low Taylor's Reynolds number Rex = 22. Without considering 
the issue of anomalous scaling in the region y < yc = 4:, they mostly concentrated on 
the cases with y > 4:, and ended up with the conclusion that for ?/ > 4 the properties of 
the statistics are not universal, but varies with the forcing spectrum slope. It is however 
difficult to consider these as conclusive results, because of the large error bars affecting 
the data. Later, in Ref. [2Z], another numerical study of the same problem at a much 
higher resolution (up to Rex = 220) has given some support, even if mostly based on 
semi-quantitative results, to the scenario drawn in the previous section in the context of 
passive transport. In particular, we report here the behaviour of longitudinal velocity 
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Figure 6. Left: PDF of the velocity increments, for y = 6.0, for three separations 
ri = 3477 and r2 — lAr] in the inertial range, and — WArj in the energy containing 
range. Distributions are normahsed to have unit variance. Right: the same as in the left 
panel, but for the case y = 3.5. The incompressible Navier-Stokes equations are solved 
with a second-order hyper-viscous dissipative term oc z/A^. Temporal integration has 
been carried over for about 20 — 30 large-eddy turnover times. DNS data refer to 256'^ 
simulations. The range of the forcing, in Fourier space, extends down to the maximum 
resolved wavenumber. 

increments PDFs for two among the different runs done in [2Z|, with y = 3.5 and y = Q. 
For y = 6 > tjc, a.s shown in Fig. IHl left panel, the usual intermittent behaviour for the 
PDFs is found: the curves P{Srv) at different separations r do not rescale one onto 
the other. On the contrary, for y = 3.5 < yc, the PDFs at various scales are almost 
indistinguishable (see Fig. El right panel): a signature of the absence of intermittent 
corrections at least for the core of the distribution. It is worth noticing that this result 
is far from trivial, because it coincides with the RG predictions EH 1221 1231 121j in the 
region where RG calculations are well beyond their range of validity. These findings, 
even if qualitative, fit rather well with those of the passive scalar case (see Fig. 0]). 

We may try to push forward this indication in the light of the theory for linear 
systems. For the 3d Navier-Stokes dynamics, the stationary equations for multi-point 
correlators Cp{r,t) = {Vai{ri)va2{r2) . . - Vapirp)) can be sketched as, 

Tp+i Cp+i + uDpCp + T2CP-2 = 0, (10) 

where Fp+i is the integro-differential linear operator coming from the inertial and 
pressure terms. Dp is the differential operator describing dissipative effects and JF2 is 
the two-point forcing correlator. Unfortunately, since the hierarchy (|TO|l is unclosed, 
no straight-forward analytical approache can be done. Even though nothing can be 
rigorously said, one may still be tempted to assume that anomalous scaling is brought 
by the inertial operator. However, the presence of finite size effects, due to the limited 
inertial range, in DNS data of 3d Navier-Stokes simulations j2Z| provide only a semi- 
quantitative support to this scenario. As we shall see in the next section, more 
quantitative statements can be made in the context of shell models for turbulence, 
for which resolution constraints are much less severe. 
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3.2. Shell models for turbulence 

Shell models of turbulence are dynamical systems mimicking Navier-Stokes nonlinear 
evolution [SE]- Their main advantage relies on the possibility of performing high 
Reynolds number simulations to measure statistical properties, say intermittent 
corrections, with high accuracy. The model we use, proposed in [SZ|, is an improved 
version of the GOY model jHHl IHl] (see also for a recent review). The evolution 
equation is : 

= i{knUl_2Un-l + h kn-lU*,_iUn+l + (1 + 6) /c„-2?^n+lMn+2) " ^klun + /„, • (H) 

The velocity field fluctuation at the wavenumber fc„, with kn = '2'^ko, is expressed in 
terms of the complex variable m„; 6 is a free parameter and u indicates the viscosity. 
The number of shells varies as n = {0, . . . , N^ax}- Remarkably, for large scale forcing, 
the structure functions show anomalous scaling: 

deviating from the dimensional prediction (p = p/3. Moreover, anomalous scaling 
exponents are found to be universal with respect to the forcing, provided it is large 
scale (see [10]). 

We aim at investigating the properties of the model stirred by a white-in-time, 
Gaussian field, with zero mean and spectrum 

{\W = fokt (13) 

being /o the forcing intensity. A similar injection mechanism was proposed in |3I] where 
authors tried to compare the power-law forced shell models with fractal-grid induced 
turbulence In order to understand the effect of forcing as (fT3|l let us start with the 
energy flux, which obeys an exact equation. Denoting with En = ^ I]^=o(l'""P/2) the 
time derivative of the energy content up to the scale A^, the energy balance relation can 
be written as : 

N N 

En + i^Y. kl{\ul\) = -k^UN + E(l/nH • (14) 

n=0 n=0 

In the previous expression, Un = (Tn), where 7/v is the energy flux through the shell 
defined as 

Tn = ^ ^^U%ulf^^UN+2 + ^ Un-iU%Un+1^ ■ (15) 

For the forcing p3|l and for u 0, we get an exact prediction in the stationary state: 

= Aik-' + Bik~'+\ (16) 

where the constants Ai,Bi, bound by the equation of motion to have opposite signs, 
depend on the forcing details (/o, (3) and on the integral wavenumber k^. At this level, 
where everything is exact, we see that when (3 approaches the critical value (/? = 0), a 
large number of shells is needed to distinguish leading terms from subleading ones in 
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Figure 7. Left: Shell fluxes measured in the presence of a power-law forcing. In 
our runs, we adopted the following choice of parameters: b — —0.4, ly — 10^^^, 
Nmax — 40. For the random forcing, acting on all shells of the inertial range, we 
chose the following values of /?: f3 — 0.3 (green squares), (3 — —0.15 (red circles). Fits 
are made with a superposition of two power-law contributions, in agreement with H16|l . 
Insets: compensation with the leading power law (dotted) or with two power laws for 
(3 — 0.3 (bottom inset) and /3 = —0.15 (top inset). Right: same as left panel but for 
the second order flux moment, to be compared with 120(1 . Also here the two power 
laws terms have opposite signs. 



the third order moment. In our runs, the choice of a large value for Nmax{= 40), allow 
us to have a good control of such effects. 

Similarly to the NS case ()10|). we can write the unclosed hierarchy for the multi- 
point correlators of shell variables: 

^Cp(n) = Mp+i(n,ri') Cp+i(n') + F2(ii,n') Cp_2(ll') , (17) 

where Cp{n) = {uniUn2---Unp) is the generic correlation function of order p, Mp+i{n, nf) is 
the linear operator coming from the inertial terms, and F2{n,n/) is the two-point forcing 
operator. The main difference with the case of passive scalar transport, as pointed out 
in the previous subsection, is that now the system (fTTj) is not closed, i.e. we have 
more unknowns than equations. Without any ambition of a rigorous approach, we may 
imagine that the general solution is characterised by two terms. The first, Cp, given by 
the dimensional matching between the inertial operator and the forcing operator, 

Mp+i(n,n') Cj+i(7i') = F2(n,n') C^.^inf). (18) 

The second, Zp, associated to a "zero mode" of the inertial operator, namely a solution 
of the homogeneous equation, Mp+i(n, n')Zp+i(n') = . Here, as for the hnear problem, 
the non homogeneous terms are the only contributions depending on the forcing, while 
the homogeneous ones, dimensionally unconstrained, may show anomalous scaling. In 
particular, for the structure functions Sp{kn), we have 

Sp{K) ~ k~^- + kn^^^"^^"^ , (19) 
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where Cp^'^il^) = p/3(l — 13). This imphes the following behaviour for the second order 
moment of the flux (equivalent to the sixth-order structure function): 

K = {r^) = A,k-'-' + B,kT^^-'\ (20) 

where the power 1.8 coincides with the anomalous exponent measured for sixth-order 
structure function with a large scale forcing, and the value 2/3 (/3 — 1) comes from the 
dimensional prediction (fT^ . 

In Fig. [7| the first and second moments of the flux variables are shown for both the 
subleading and dominant forcing with (3 = —0.15 and /? = 0.3, respectively. It is easy 
to see that to get a good agreement of the numerical results with the prediction (pUj) . 
both the leading and the subleading power laws have to be taken into account. For both 
cases with P = 0.3 and /3 = —0.15, at least for the order of the statistics we could reach, 
the anomalous scaling exponents have values in agreement with those found for a large 
scale forcing, supporting the universality scenario. 

4. Conclusions 

We have discussed the problem of small scale fluctuations in turbulent systems stirred 
at all scales by a power-law forcing. The main question in our study concerns anomalous 
scaling and universality of scaling exponents. 

For linear systems - i.e. passive transport-, we find clear evidence for the 
universality of anomalous fluctuations and our results fit well in the zero mode 
scenario. In the regime in which forcing is subleading anomalous scaling is recovered in 
quantitative agreement with the case of large scale forcing. In the forcing dominated 
regime, the dimensional scaling imposed by the injection mechanism overwhelms the 
anomalous fluctuations of low order moments. Nevertheless, anomalous scaling shows 
up again for high enough moments, independently of the forcing spectrum slope, as 
confirmed by the existence of a unique saturation exponent in all regimes. 

The results obtained for nonlinear turbulent systems point in the same direction. 
Indeed both the semi-quantitative results obtained in the context of 3d Navier-Stokes 
turbulence and those more quantitative issuing from the investigation of the shell models 
are compatible with the linear systems scenario for universality. However, in the presence 
of power-law forcing acting in the turbulent energy cascade range, as shown both in the 
linear and nonlinear cases, scaling properties are strongly spoiled by the beating of 
leading and subleading terms. This effect is particularly strong due to cancellations 
induced by different signed pre-factors in the power-law terms. 

Strong cancellation effects, which apparently are a mere technical question, may 
lead to misinterpretation in analysing data. We wonder, for example, if the observed 
multi-fractal behaviour in the one- dimensional Burgers equation stirred by a power-law 
forcing might be a spurious effect ji^ . 
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